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Starting from the general Horndeski action, we derive the most general effective theory for scalar 
perturbations around flat space that allows to screen fifth forces via the Vainshtein mechanism. The 
effective theory is described by a generalization of the Galileon Lagrangian, which we use to study the 
stability of spherically symmetric configurations exhibiting the Vainshtein effect. In particular, we 
discuss the phenomenological consequences of a scalar-tensor coupling that is absent in the standard 
Galileon lagrangian. This coupling controls the superluminality and stability of fluctuations inside 
the Vainshtein radius in a way that depends on the density profile of a matter source. Particularly 
we find that the vacuum solution is unstable due to this coupling. 
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Introduction: The late time acceleration of the Uni- 
verse has become one of the most challenging problems 
in theoretical physics. The cosmological constant is the 
simplest option but its required value is unacceptably 
small from the point of view of effective field theory. Var- 
ious models have been proposed as alternatives to a pure 
cosmological constant. Many of them make use of scalar 
fields to achieve the accelerated expansion. Scalar de- 
grees of freedom also appear in modified gravity models 
where the acceleration is attributed to the failure of Gen- 
eral Relativity (GR) on cosmological scales. However, in 
general scalars mediate fifth forces strongly constrained 
by precision tests of gravity at solar system scales: any 
dark energy and modified gravity model involving scalar 
fields should accommodate a mechanism to suppress the 
scalar interaction on small scales, where stringent con- 
straints on the deviation from GR are applicable. 

There has been interesting progress in developing 
screening mechanisms to suppress scalar interactions on 
small scales. One of the oldest ideas is the Vainshtein 
mechanism [1], originally discovered in the context of 
massive gravity (see Ref. [2] for a review). Massive gravi- 
tons have five polarisations instead of the two of GR: the 
helicity-0 mode couples to matter and modifies the pre- 
diction of GR. In a linear approximation, this helicity-0 
mode does not decouple in the massless limit, leading 
to the so-called vDVZ discontinuity [3, 4]. This prob- 
lem can be solved by the Vainshtein mechanism. If the 
graviton mass is small, the derivative self-interactions of 
the helicity-0 mode become important at much larger 
distances compared with the Schwarzschild radius of a 
source, and they suppress the coupling of the helicity-0 
mode to matter. The essential feature of this mechanism, 
which makes it different from other proposals, is that 
derivative self-interactions of the scalars are the main re- 
sponsible for the screening mechanism, without requiring 



any particular form of scalar potential, and any couplings 
of the scalar to matter other than gravitational. In the 
more direct extension of the Fierz-Pauli massive gravity 
[5] coupled to a source, the non-linear interactions in- 
volve higher derivatives in the equations of motion, and 
this leads to the Boulware-Deser ghost [6]. The non- 
linear interactions that give a scalar equation of motion 
containing only up to second derivatives have been con- 
structed imposing a Galilean symmetry, that allows only 
a handful of terms in the scalar lagrangian [7]. The co- 
variant version of this theory has been obtained as the 
most general scalar tensor theory whose equation of mo- 
tion contains up to the second derivative of the fields 
[8, 9]. The action describing this theory was indeed first 
obtained by Horndeski in 1974 [10] and thus is now called 
the Horndeski action. 

In this letter, we discuss the general effective theory 
that describes the Vainshtein mechanism in flat space- 
time starting from the Horndeski action. Although we 
derive this effective theory from it, the infrared comple- 
tion of the theory can be described by a different ac- 
tion. Indeed the effective theory we analyze includes 
suitable decoupling limits of other scenarios, such as 
Dvali-Gabadadze-Porrati (DGP) braneworld model [11] 
or ghost-free massive gravity [12]. 

The Horndeski action: Our starting point is the Horn- 
deski action [8, 10], given by 
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The first term controls the dynamics of the scalar field, 
and contains the following four pieces: 
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with X := —g fxu d fl (j)d l ,(t)/2, and (V M V„</>) n a short nota- 
tion for trace of the n-time contraction of V^V^^. The 
notation As ■= dA/dB will be used throughout the text. 
Notice that, in the second part of the action (1), we are 
choosing a frame in which matter does not directly couple 
to the scalar, but only to gravity. 

An effective theory for the Vainshtein mechanism: 

For the purpose of studying the Vainshtein mechanism, 
we consider a Minkowski background with a constant 
scalar field <fi = 4>q and study the deviations around it, 
namely 



^0 + 7T, = Tlfa, + 
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The dynamics of fluctuations around this background 
controls the behavior of gravitational interactions. No- 
tice that for the background to be a solution, one needs 
K(4>o) = K $((/>, 0) = 0. We construct an effective theory 
for the dynamics of fluctuations aimed to describe the 
Vainshtein mechanism. As explained above, the distinc- 
tive feature of the Vainshtein mechanism is that deriva- 
tive self-interactions of the scalar are sufficient to con- 
trol its behavior around a source, and screen its effect. 
Therefore, in order to identify these derivative interac- 
tions and realize the Vainshtein mechanism as described 
above, we demand that the scalar does not directly cou- 
ple to the source, leaving possible couplings to matter 
through mixing with gravity. In order to do so, we im- 
pose that the scalar action respects the Galileon symme- 
try tt —t tt + c + b^x^, introduced in [7]. This symmetry 
prevents direct couplings between the tt fluctuation and 
the energy momentum tensor. Moreover, the symmetry 
does not allow a potential for the scalar field tt, and cou- 
plings between first derivatives of tt and tt itself. As a 
consequence, it is the symmetry that we need in order 
to isolate the pieces of the action that depend purely on 
derivative self-interactions of the scalar sector. 

We achieve this effective action by expanding the Horn- 
deski action (1) in terms of the fluctuations (7), using the 
following assumptions: the fields tt and are small, 
hence we neglect higher order interactions containing the 
tensor perturbations h^u, as well as terms containing 
higher order powers of the scalar fluctuation tt and its 
first derivatives only. On the other hand, we keep all 
terms with second order derivatives of tt that preserve 
the Galileon symmetry, and will provide the necessary 
self-interactions to realize the Vainshtein screening mech- 
anism. In other words, this means to Taylor-expand the 
functions K, and Gi to first order in the quantities X and 
<j>, and the curvature tensors G^ v and R to first and sec- 
ond order in ha V , respectively: all the higher derivatives 
of these functions are required to be suppressed, in order 



to obtain an effective action that preserves the Galilean 
symmetry. By doing this expansion, we schematically 
find the following terms 
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The kinetic term for metric perturbations is contained 
in Gi{4>, x)R. Thus we introduce the Planck scale as 
G4 = Mpi/2 and canonically normalise as h^u = 
Mpihfiu. Wc further introduce a new mass dimension A 
so that terms involving second order derivatives of tt have 
the appropriate dimensions. By explicitly expanding the 
Horndeski action (I) under these prescriptions, we find 
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where £ gal are the so-called Galileon terms [7] 
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and the tensor-scalar couplings Xi are given by 
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where the short hand notations n M „ = d^d v TT, II" = 
n MQ n Q ^ ...n A „ and [IP] = II" ^ are introduced and 
£ a P is the linearised Einstein tensor. We defined 
seven dimcnsionlcss parameters £, 77, fi, v, w, a, ft as fol- 
lows [13] 
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where all functions are evaluated at the background, 
(f> = 4>o and X = 0. We assume these dimensionless 
parameters are 0(1) when they are non-zero. 

Notice that, as anticipated, the above action respects 
the Galileon symmetry (up to total derivative terms). 
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On the other hand, the coupling between the metric and 
scalar perturbations h^xj^J and h^ u xj?J can be elimi- 
nated by the local field redefinition [15] 

- 2a 

hp* = hp U - 2£7T77 M „ + -j^dfj,Ttd v 'K, (13) 



so that the action (9) now becomes 
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However, the coupling W 'Xfaj cannot be removed by a 
local field redefinition [15]. The transformation (13) in- 
troduces a coupling between 7r and T„ v via the de-mixing 
of scalar from gravity. 

This is the most general form of the effective theory 
aimed to analyse the Vainshtein mechanism, which in- 
volves up to the second derivatives of the fields in the 
equations of motion. We should emphasise that, in our 
discussion, the scalar field perturbation 7r and its first 
derivative 8tt are assumed to be small. There are alter- 
native screening mechanisms using the non-linearity of 



the field itself or the first derivative (see Rcf. [14] for a 
review). In these cases, we should expand the functions 
K and Gi at all orders, hence it is required to know the 
functional forms of these functions. This information is 
not needed in our approach, which isolates the derivative 
self-interactions that are needed to identify and analyze 
the Vainshtein mechanism in its traditional form. 

As wc mentioned before, although we reached this re- 
sult from the Horndcski action, the infrared completion 
of the theory can be different from the Horndeski ac- 
tion. For example, the recently proposed massive grav- 
ity theory reduces to (14) in the decoupling limit with 
■q = fj, = u = m = [15]. In this case, A = m 2 Mp\, with 
m being the graviton's mass, while a and /3 are related to 
the two dimensionless parameter of this massive gravity 
theory, which are usually referred as and a^. A further 
example is the decoupling limit of the DGP braneworld 
model, which is reproduced hy a = ft = v = m = [16]. 



In this case, A 



Mpi where r c is the cross-over scale 



between four-dimensional and five-dimensional gravity. 
In both models, tt is not a fundamental scalar field, but 
it is the helicity-0 component of a graviton. 

Spherically symmetric solutions: From the effective 
action (14), it is easy to find static spherically symmetric 
solutions. We consider the following configuration 

ha = -2$, kj = -2*%, = + 7r(r). (15) 
The field equations yield 



P(x, A) ee £A{r) + (| + 3£ 2 ) x + (/z + 6a£ - 3/3A(r)) x 2 + (v + 2a 2 + 4/3£) x 3 - 3f3 2 x 5 = 0, (16) 
y (r) = I3x 3 +A{r) 1 (17) 



where wc define 
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The function M(r) represents the mass of a spherically 
symmetric, pressureless matter source, up to a radius r. 
Outside the surface r s of the matter source, A(r) can be 
written as 
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where ry is the Vainshtein radius depending on the total 
mass M of the source. Thus, well inside the Vainshtein 
radius, A(r) 3> 1 and the non-linearity of x becomes 
large. For the Sun, the Vainshtein radius ry is O.lkpc 
if A is associated with the energy scale for the current 



accelerated expansion of the Universe A ~ Mp\H 2 . So- 
lutions to these equations were studied in [17-19] within 
massive gravity models, and the analysis was extended 
to the general case in [13, 20]. 

Stability of these solutions can be studied by consid- 
ering small perturbations around the background it = 
7r (r) + tp(t,r,Q,),$ = $(r) + <S$(i,r, fi), * = *(r) + 
8&(t,r,n) where Q represents angular coordinates. By 
expanding the effective action up to the second order in 
these small perturbations we find 
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In addition to the quadratic action for scalar perturba- 
tions, there is also a coupling between metric pertur- 
bations and the scalar perturbations. The stability of 
fluctuations depend strongly on the presence of two par- 
ticular couplings in the Lagrangian (14), one given by 
9 Al 7r9 !/ 7rT' iI/ and the other by h^ v Xj^J , which have a 
and (3 as coefficients, respectively. Some of their effects 
have already been analysed in the literature, but other 
have not. We will consider their consequences in what 
follows, by analyzing each relevant case separately. 

Case 1: a = (5 = 0. In this case the action (14) re- 
duces to the Galileon action introduced by Ref. [7]. The 
stability conditions were studied in detail in the same ref- 
erence, hence we do not repeat them here. Interestingly, 
it was found that the stability conditions force the the 
speed of sound for fluctuations propagating radially to 
be greater than one. On the other hand, if the fourth 
order Galileon term Cf 11 is included, then Kq <C K t ,K r 
and the propagation of angular fluctuations is extremely 
sub-luminal, invalidating a quasi-static approximation. 

Case II: a ^ and (3—0. A new feature of this class of 
models is the disformal coupling between matter and the 
scalar, namely d^nd^TrT^ . As pointed out in Ref. [21], 
this coupling has deep implications for the stability. In- 
side a matter source, the coupling introduces a kinetic 
term proportional to ap(d t f) 2 ■ If a < 0, the fluctuations 
behave as ghosts, forcing us to choose a > 0. In the case 
of massive gravity, there is only one dimensionless param- 
eter, a = l + 3a3 (in the theory where /3 = a^+Aa^ = 0), 
leading to the so-called restricted Galileon [21]. In this 
case, there is no Vainshtein solution that can be con- 
nected to the asymptotically flat solution. Around the 
non-flat asymptotic solution, the radial propagation is 
sub-luminal, even though the extreme sub-luminarity of 
the angular propagation persists [21]. 

Case III: a ^ 0, /3 ^ 0. In this case, the coupling be- 
tween the metric and the scalar WX^J gives a different 
picture from the previous cases and this has not been 
analysed so far. Inside the Vainshtein radius the solution 
that reduces to GR near a matter source only exists if 
j3 > 0, and the solution is given by [18, 19] 



.T = ±1 
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y ~ A(r) 3> 1 contribution in the limit r <C ry, as shown 
from its definition (18). In the limit A(r) ^> 1, we can 
also ignore the coupling between the scalar and metric 
perturbations. By taking into account only the contribu- 
tions depending on A = (ry/r) 3 outside the surface of 
the source r > r s , we find in this limit 

K t = 0, K r = -12pA(r)x , K n = 6pA(r)x . (26) 

Thus inside the Vainshtein radius r <C ry, but outside 
the source surface, the speed of the fluctuations are al- 
ways supcrluminal, for both radial and angular direc- 
tions. Moreover, given the fact that Kq and K r have 
opposite sign, all solutions are unstable. 

On the other hand, the instability can be avoided inside 
the source. For definiteness, we consider a distribution 
of matter given by p(r) = po(r/r s ) s for r < r s . Inside 
the matter source r < r s , the Ki functions are given by 

K t = 4(a + 3(3x )(3 + s)A(r), 

K r = -12(3x A(r), K n - -6(2 + s)f3x A(r), (27) 

where A(r) = p(r)/ [2(3 + s)M P iA 3 ] . The stability and 
superluminality depend on the slope of the density pro- 
file, s, as well as on the values of the coefficients a, f3. 
The speeds of sound for radial and angular fluctuations 
read 



12f3x 



4(a + 3/3x )(3 + ,s) : 

6(2 + s)(3x 
4(a + 3/3x )(3 + s) : 



(28) 
(29) 



Since x is constant and A(r) ^> 1 inside the Vainshtein 
radius, all functions Ki in (22) are dominated by the 



and their ratio is c 2 /^ = 2/(2 + s). Choosing the neg- 
ative branch in eq. (25), s > —2, and a positive and 
sufficiently large, both the previous expressions can be 
made positive avoiding instabilities. For appropriate val- 
ues of a and s, moreover, the perturbations can be made 
subluminal, but not so subluminal to invalidate a quasi- 
static approximation. 

Discussion: In this letter, we derived the effective ac- 
tion that describes the Vainshtein mechanism from the 
Horndcski action. Using the effective theory, we stud- 
ied spherically symmetric solutions and perturbations 
around them. The solutions can be classified by two pa- 
rameters a and (3. If a ^ 0, there appears a disformal 
coupling between the scalar and energy-momentum ten- 
sor. If (3 7^ 0, the coupling between metric and scalar 
cannot be diagonalised by a local field transformation. 
We found that if (3 ^ 0, the vacuum spherically symmet- 
ric solution is pathological. Inside the Vainshtein radius, 
at the leading order in (ry /r) 3> 1, the temporal kinetic 
term vanishes, leading to an extreme superluminal prop- 
agation. Moreover, the spatial kinetic term for radial 
and angular modes have an opposite sign, so the pertur- 
bations are unstable. This instability can be cured inside 
a matter source, choosing an appropriate density profile. 
In order to avoid the instability of vacuum solutions, on 
the other hand, we would need to fine-tune (3 such that 
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p < 0((?' s /ry) 2 ). For the case of the Earth, this condi- 
tion would imply j3 < 1CP 22 for a cut-off scale of order 
A ~ Mp\Hq. This imposes a strong restriction on the 
function G5 appearing in the Hordenski action. 

For a further understanding of this coupling mediated 
by /3 7^ 0, it would be interesting to solve the system nu- 
merically when the spherical symmetry is relaxed, using 
the techniques developed in Rcf. [22]. Due to the com- 
plex tensor structure of the equations, screening works 
very differently in non-spherically symmetric system [23] . 
Another interesting possibility is to consider the effect 
of pressure, since it has already been pointed out that 
pressure gives contributions to K r and Kq through the 
disformal coupling depending on a [21, 24]. 

Recently, it was argued that, except for the model with 
only the cubic Galileon term £| a , the theory looses its 
predictive power due to quantum corrections at C(l)mm 
if A is associated with the energy scale for the current 
accelerated expansion of the Universe A ~ Mp\Hq [25]. 
There is still a debate whether the disformal coupling is 
enough to cure this problem [21, 26] but we should point 



out that the situation might be worse when we consider a 
spherically symmetric vacuum solution as a background, 
if /3 7^ 0. In this case, the temporal kinetic term is order 
unity K t ~ (9(1), which is much smaller than the one 
reduced from an order estimation K t ~ (ry/r) 3 . Again 
the existence of matter significantly increases the kinetic 
term K t ~ (p/Mp\A 3 ) and the strong coupling scale de- 
pends on the matter distribution. 
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